Risk of rapid evolutionary escape from biomedical interventions targeting SARS-CoV-2 spike protein.

The spike protein receptor-binding domain (RBD) of SARS-CoV-2 is the molecular target for many vaccines and antibody-based prophylactics aimed at bringing COVID-19 under control. Such a narrow molecular focus raises the specter of viral immune evasion as a potential failure mode for these biomedical interventions. With the emergence of new strains of SARS-CoV-2 with altered transmissibility and immune evasion potential, a critical question is this: how easily can the virus escape neutralizing antibodies (nAbs) targeting the spike RBD? To answer this question, we combined an analysis of the RBD structure-function with an evolutionary modeling framework. Our structure-function analysis revealed that epitopes for RBD-targeting nAbs overlap one another substantially and can be evaded by escape mutants with ACE2 affinities comparable to the wild type, that are observed in sequence surveillance data and infect cells in vitro. This suggests that the fitness cost of nAb-evading mutations is low. We then used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to vaccines or other prophylactics that rely on one or two antibodies for protection can develop quickly -and repeatedly- under positive selection. Predicted resistance timelines are comparable to those of the decay kinetics of nAbs raised against vaccinal or natural antigens, raising a second potential mechanism for loss of immunity in the population. Strategies for viral elimination should therefore be diversified across molecular targets and therapeutic modalities.

SARS-CoV-2 Membrane Protein: From Genomic Data to Structural New Insights.

Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) is composed of four structural proteins and several accessory non-structural proteins. SARS-CoV-2's most abundant structural protein, Membrane (M) protein, has a pivotal role both during viral infection cycle and host interferon antagonism. This is a highly conserved viral protein, thus an interesting and suitable target for drug discovery. In this paper, we explain the structural nature of M protein homodimer. To do so, we developed and applied a detailed and robust in silico workflow to predict M protein dimeric structure, membrane orientation, and interface characterization. Single Nucleotide Polymorphisms (SNPs) in M protein were retrieved from over 1.2 M SARS-CoV-2 genomes and proteins from the Global Initiative on Sharing All Influenza Data (GISAID) database, 91 of which were located at the predicted dimer interface. Among those, we identified SNPs in Variants of Concern (VOC) and Variants of Interest (VOI). Binding free energy differences were evaluated for dimer interfacial SNPs to infer mutant protein stabilities. A few high-prevalent mutated residues were found to be especially relevant in VOC and VOI. This realization may be a game-changer to structure-driven formulation of new therapeutics for SARS-CoV-2.

Nebulized delivery of a broadly neutralizing SARS-CoV-2 RBD-specific nanobody prevents clinical, virological, and pathological disease in a Syrian hamster model of COVID-19.

There remains an unmet need for globally deployable, low-cost therapeutics for the ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Previously, we reported on the isolation and in vitro characterization of a potent single-domain nanobody, NIH-CoVnb-112, specific for the receptor-binding domain (RBD) of SARS-CoV-2. Here, we report on the molecular basis for the observed broad in vitro neutralization capability of NIH-CoVnb-112 against variant SARS-CoV-2 pseudoviruses. The structure of NIH-CoVnb-112 bound to SARS-CoV-2 RBD reveals a large contact surface area overlapping the angiotensin converting enzyme 2 (ACE2) binding site, which is largely unencumbered by the common RBD mutations. In an in vivo pilot study, we demonstrate effective reductions in weight loss, viral burden, and lung pathology in a Syrian hamster model of COVID-19 following nebulized delivery of NIH-CoVnb-112. These findings support the further development of NIH-CoVnb-112 as a potential adjunct preventative therapeutic for the treatment of SARS-CoV-2 infection.Abbreviations: ACE2 - angiotensin converting enzyme 2BSA - buried surface areaCDR - complementary determining regionRBD - receptor binding domainRBM - receptor-binding motifSARS-CoV-2 - severe acute respiratory syndrome coronavirus 2.

Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2.

The most recent Omicron variant of SARS-CoV-2 has caused global concern and anxiety. The only thing certain about this strain, with a large number of mutations in the spike protein, is that it spreads quickly, seems to evade immune defense, and mitigates the benefits of existing vaccines. Based on the ultra-large-scale ab initio computational modeling of the receptor binding motif (RBM) and the human angiotensin-converting enzyme-2 (ACE2) interface, we provide the details of the effect of Omicron mutations at the fundamental atomic scale level. In-depth analysis anchored in the novel concept of amino acid-amino acid bond pair units (AABPU) indicates that mutations in the Omicron variant are connected with (i) significant changes in the shape and structure of AABPU components, together with (ii) significant increase in the positive partial charge, which facilitates the interaction with ACE2. We have identified changes in bonding due to mutations in the RBM. The calculated bond order, based on AABPU, reveals that the Omicron mutations increase the binding strength of RBM to ACE2. Our findings correlate with and are instrumental to explain the current observations and can contribute to the prediction of next potential new variant of concern.

Evolution of the SARS-CoV-2 spike protein in the human host.

Recently emerged variants of SARS-CoV-2 contain in their surface spike glycoproteins multiple substitutions associated with increased transmission and resistance to neutralising antibodies. We have examined the structure and receptor binding properties of spike proteins from the B.1.1.7 (Alpha) and B.1.351 (Beta) variants to better understand the evolution of the virus in humans. Spikes of both variants have the same mutation, N501Y, in the receptor-binding domains. This substitution confers tighter ACE2 binding, dependent on the common earlier substitution, D614G. Each variant spike has acquired other key changes in structure that likely impact virus pathogenesis. The spike from the Alpha variant is more stable against disruption upon binding ACE2 receptor than all other spikes studied. This feature is linked to the acquisition of a more basic substitution at the S1-S2 furin site (also observed for the variants of concern Delta, Kappa, and Omicron) which allows for near-complete cleavage. In the Beta variant spike, the presence of a new substitution, K417N (also observed in the Omicron variant), in combination with the D614G, stabilises a more open spike trimer, a conformation required for receptor binding. Our observations suggest ways these viruses have evolved to achieve greater transmissibility in humans.

Virus-specific editing identification approach reveals the landscape of A-to-I editing and its impacts on SARS-CoV-2 characteristics and evolution.

Upon SARS-CoV-2 infection, viral intermediates specifically activate the IFN response through MDA5-mediated sensing and accordingly induce ADAR1 p150 expression, which might lead to viral A-to-I RNA editing. Here, we developed an RNA virus-specific editing identification pipeline, surveyed 7622 RNA-seq data from diverse types of samples infected with SARS-CoV-2, and constructed an atlas of A-to-I RNA editing sites in SARS-CoV-2. We found that A-to-I editing was dynamically regulated, varied between tissue and cell types, and was correlated with the intensity of innate immune response. On average, 91 editing events were deposited at viral dsRNA intermediates per sample. Moreover, editing hotspots were observed, including recoding sites in the spike gene that affect viral infectivity and antigenicity. Finally, we provided evidence that RNA editing accelerated SARS-CoV-2 evolution in humans during the epidemic. Our study highlights the ability of SARS-CoV-2 to hijack components of the host antiviral machinery to edit its genome and fuel its evolution, and also provides a framework and resource for studying viral RNA editing.

A Review of Alpha-1 Antitrypsin Binding Partners for Immune Regulation and Potential Therapeutic Application.

Alpha-1 antitrypsin (AAT) is the canonical serine protease inhibitor of neutrophil-derived proteases and can modulate innate immune mechanisms through its anti-inflammatory activities mediated by a broad spectrum of protein, cytokine, and cell surface interactions. AAT contains a reactive methionine residue that is critical for its protease-specific binding capacity, whereby AAT entraps the protease on cleavage of its reactive centre loop, neutralises its activity by key changes in its tertiary structure, and permits removal of the AAT-protease complex from the circulation. Recently, however, the immunomodulatory role of AAT has come increasingly to the fore with several prominent studies focused on lipid or protein-protein interactions that are predominantly mediated through electrostatic, glycan, or hydrophobic potential binding sites. The aim of this review was to investigate the spectrum of AAT molecular interactions, with newer studies supporting a potential therapeutic paradigm for AAT augmentation therapy in disorders in which a chronic immune response is strongly linked.

Ultrapotent neutralizing antibodies against SARS-CoV-2 with a high degree of mutation resistance.

Many SARS-CoV-2 neutralizing antibodies (nAbs) lose potency against variants of concern. In this study, we developed 2 strategies to produce mutation-resistant antibodies. First, a yeast library expressing mutant receptor binding domains (RBDs) of the spike protein was utilized to screen for potent nAbs that are least susceptible to viral escape. Among the candidate antibodies, P5-22 displayed ultrahigh potency for virus neutralization as well as an outstanding mutation resistance profile. Additionally, P14-44 and P15-16 were recognized as mutation-resistant antibodies with broad betacoronavirus neutralization properties. P15-16 has only 1 binding hotspot, which is K378 in the RBD of SARS-CoV-2. The crystal structure of the P5-22, P14-44, and RBD ternary complex clarified the unique mechanisms that underlie the excellent mutation resistance profiles of these antibodies. Secondly, polymeric IgG enhanced antibody avidity by eliminating P5-22's only hotspot, residue F486 in the RBD, thereby potently blocking cell entry by mutant viruses. Structural and functional analyses of antibodies screened using both potency assays and the yeast RBD library revealed rare, ultrapotent, mutation-resistant nAbs against SARS-CoV-2.

The stalk domain of SARS-CoV-2 NSP13 is essential for its helicase activity.

COVID-19, caused by SARS-CoV-2, has been spreading worldwide for more than two years and has led to immense challenges to human health. Despite the great efforts that have been made, our understanding of SARS-CoV-2 is still limited. The viral helicase, NSP13 is an important enzyme involved in SARS-CoV-2 replication and transcription. Here we highlight the important role of the stalk domain in the enzymatic activity of NSP13. Without the stalk domain, NSP13 loses its dsRNA unwinding ability due to the lack of ATPase activity. The stalk domain of NSP13 also provides a rigid connection between the ZBD and helicase domain. We found that the tight connection between the stalk and helicase is necessary for NSP13-mediated dsRNA unwinding. When a short flexible linker was inserted between the stalk and helicase domains, the helicase activity of NSP13 was impaired, although its ATPase activity remained intact. Further study demonstrated that linker insertion between the stalk and helicase domains attenuated the RNA binding ability and affected the thermal stability of NSP13. In summary, our results suggest the crucial role of the stalk domain in NSP13 enzymatic activity and provide mechanistic insight into dsRNA unwinding by SARS-CoV-2 NSP13.

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN.

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

An extended motif in the SARS-CoV-2 spike modulates binding and release of host coatomer in retrograde trafficking.

ß-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for αCOPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. αCOPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking.

The low abundance of CpG in the SARS-CoV-2 genome is not an evolutionarily signature of ZAP.

The zinc finger antiviral protein (ZAP) is known to restrict viral replication by binding to the CpG rich regions of viral RNA, and subsequently inducing viral RNA degradation. This enzyme has recently been shown to be capable of restricting SARS-CoV-2. These data have led to the hypothesis that the low abundance of CpG in the SARS-CoV-2 genome is due to an evolutionary pressure exerted by the host ZAP. To investigate this hypothesis, we performed a detailed analysis of many coronavirus sequences and ZAP RNA binding preference data. Our analyses showed neither evidence for an evolutionary pressure acting specifically on CpG dinucleotides, nor a link between the activity of ZAP and the low CpG abundance of the SARS-CoV-2 genome.

Synthetic virions reveal fatty acid-coupled adaptive immunogenicity of SARS-CoV-2 spike glycoprotein.

SARS-CoV-2 infection is a major global public health concern with incompletely understood pathogenesis. The SARS-CoV-2 spike (S) glycoprotein comprises a highly conserved free fatty acid binding pocket (FABP) with unknown function and evolutionary selection advantage1,2. Deciphering FABP impact on COVID-19 progression is challenged by the heterogenous nature and large molecular variability of live virus. Here we create synthetic minimal virions (MiniVs) of wild-type and mutant SARS-CoV-2 with precise molecular composition and programmable complexity by bottom-up assembly. MiniV-based systematic assessment of S free fatty acid (FFA) binding reveals that FABP functions as an allosteric regulatory site enabling adaptation of SARS-CoV-2 immunogenicity to inflammation states via binding of pro-inflammatory FFAs. This is achieved by regulation of the S open-to-close equilibrium and the exposure of both, the receptor binding domain (RBD) and the SARS-CoV-2 RGD motif that is responsible for integrin co-receptor engagement. We find that the FDA-approved drugs vitamin K and dexamethasone modulate S-based cell binding in an FABP-like manner. In inflammatory FFA environments, neutralizing immunoglobulins from human convalescent COVID-19 donors lose neutralization activity. Empowered by our MiniV technology, we suggest a conserved mechanism by which SARS-CoV-2 dynamically couples its immunogenicity to the host immune response.

Crystal Structures of Bat and Human Coronavirus ORF8 Protein Ig-Like Domain Provide Insights Into the Diversity of Immune Responses.

ORF8 is a viral immunoglobulin-like (Ig-like) domain protein encoded by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA genome. It tends to evolve rapidly and interfere with immune responses. However, the structural characteristics of various coronavirus ORF8 proteins and their subsequent effects on biological functions remain unclear. Herein, we determined the crystal structures of SARS-CoV-2 ORF8 (S84) (one of the epidemic isoforms) and the bat coronavirus RaTG13 ORF8 variant at 1.62 Å and 1.76 Å resolution, respectively. Comparison of these ORF8 proteins demonstrates that the 62-77 residues in Ig-like domain of coronavirus ORF8 adopt different conformations. Combined with mutagenesis assays, the residue Cys20 of ORF8 is responsible for forming the covalent disulfide-linked dimer in crystal packing and in vitro biochemical conditions. Furthermore, immune cell-binding assays indicate that various ORF8 (SARS-CoV-2 ORF8 (L84), ORF8 (S84), and RaTG13 ORF8) proteins have different interaction capabilities with human CD14+ monocytes in human peripheral blood. These results provide new insights into the specific characteristics of various coronavirus ORF8 and suggest that ORF8 variants may influence disease-related immune responses.

Highly conserved s2m element of SARS-CoV-2 dimerizes via a kissing complex and interacts with host miRNA-1307-3p.

The ongoing COVID-19 pandemic highlights the necessity for a more fundamental understanding of the coronavirus life cycle. The causative agent of the disease, SARS-CoV-2, is being studied extensively from a structural standpoint in order to gain insight into key molecular mechanisms required for its survival. Contained within the untranslated regions of the SARS-CoV-2 genome are various conserved stem-loop elements that are believed to function in RNA replication, viral protein translation, and discontinuous transcription. While the majority of these regions are variable in sequence, a 41-nucleotide s2m element within the genome 3' untranslated region is highly conserved among coronaviruses and three other viral families. In this study, we demonstrate that the SARS-CoV-2 s2m element dimerizes by forming an intermediate homodimeric kissing complex structure that is subsequently converted to a thermodynamically stable duplex conformation. This process is aided by the viral nucleocapsid protein, potentially indicating a role in mediating genome dimerization. Furthermore, we demonstrate that the s2m element interacts with multiple copies of host cellular microRNA (miRNA) 1307-3p. Taken together, our results highlight the potential significance of the dimer structures formed by the s2m element in key biological processes and implicate the motif as a possible therapeutic drug target for COVID-19 and other coronavirus-related diseases.

E484K and N501Y SARS-CoV 2 spike mutants Increase ACE2 recognition but reduce affinity for neutralizing antibody.

SARS-CoV2 mutants B.1.1.7, B.1.351, and P.1 contain a key mutation N501Y. B.1.135 and P.1 lineages have another mutation, E484K. Here, we decode the effect of these two mutations on the host receptor, ACE2, and neutralizing antibody (B38) recognition. The N501Y RBD mutant binds to ACE2 with higher affinity due to improved π-π stacking and π-cation interactions. The higher binding affinity of the E484K mutant is caused due to the formation of additional hydrogen bond and salt-bridge interactions with ACE2. Both the mutants bind to the B38 antibody with reduced affinity due to the loss of several hydrogen-bonding interactions. The insights obtained from the study are crucial to interpret the increased transmissibility and reduced neutralization efficacy of rapidly emerging SARS-CoV2 VOCs.

A Routine Sanger Sequencing Target Specific Mutation Assay for SARS-CoV-2 Variants of Concern and Interest.

As SARS-CoV-2 continues to spread among human populations, genetic changes occur and accumulate in the circulating virus. Some of these genetic changes have caused amino acid mutations, including deletions, which may have a potential impact on critical SARS-CoV-2 countermeasures, including vaccines, therapeutics, and diagnostics. Considerable efforts have been made to categorize the amino acid mutations of the angiotensin-converting enzyme 2 (ACE2) receptor binding domain (RBD) of the spike (S) protein, along with certain mutations in other regions within the S protein as specific variants, in an attempt to study the relationship between these mutations and the biological behavior of the virus. However, the currently used whole genome sequencing surveillance technologies can test only a small fraction of the positive specimens with high viral loads and often generate uncertainties in nucleic acid sequencing that needs additional verification for precision determination of mutations. This article introduces a generic protocol to routinely sequence a 437-bp nested RT-PCR cDNA amplicon of the ACE2 RBD and a 490-bp nested RT-PCR cDNA amplicon of the N-terminal domain (NTD) of the S gene for detection of the amino acid mutations needed for accurate determination of all variants of concern and variants of interest according to the definitions published by the U.S. Centers for Disease Control and Prevention. This protocol was able to amplify both nucleic acid targets into cDNA amplicons to be used as templates for Sanger sequencing on all 16 clinical specimens that were positive for SARS-CoV-2.

Low-dose in vivo protection and neutralization across SARS-CoV-2 variants by monoclonal antibody combinations.

Prevention of viral escape and increased coverage against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern require therapeutic monoclonal antibodies (mAbs) targeting multiple sites of vulnerability on the coronavirus spike glycoprotein. Here we identify several potent neutralizing antibodies directed against either the N-terminal domain (NTD) or the receptor-binding domain (RBD) of the spike protein. Administered in combinations, these mAbs provided low-dose protection against SARS-CoV-2 infection in the K18-human angiotensin-converting enzyme 2 mouse model, using both neutralization and Fc effector antibody functions. The RBD mAb WRAIR-2125, which targets residue F486 through a unique heavy-chain and light-chain pairing, demonstrated potent neutralizing activity against all major SARS-CoV-2 variants of concern. In combination with NTD and other RBD mAbs, WRAIR-2125 also prevented viral escape. These data demonstrate that NTD/RBD mAb combinations confer potent protection, likely leveraging complementary mechanisms of viral inactivation and clearance.

Glycosylation and Serological Reactivity of an Expression-enhanced SARS-CoV-2 Viral Spike Mimetic.

Extensive glycosylation of viral glycoproteins is a key feature of the antigenic surface of viruses and yet glycan processing can also be influenced by the manner of their recombinant production. The low yields of the soluble form of the trimeric spike (S) glycoprotein from SARS-CoV-2 has prompted advances in protein engineering that have greatly enhanced the stability and yields of the glycoprotein. The latest expression-enhanced version of the spike incorporates six proline substitutions to stabilize the prefusion conformation (termed SARS-CoV-2 S HexaPro). Although the substitutions greatly enhanced expression whilst not compromising protein structure, the influence of these substitutions on glycan processing has not been explored. Here, we show that the site-specific N-linked glycosylation of the expression-enhanced HexaPro resembles that of an earlier version containing two proline substitutions (2P), and that both capture features of native viral glycosylation. However, there are site-specific differences in glycosylation of HexaPro when compared to 2P. Despite these discrepancies, analysis of the serological reactivity of clinical samples from infected individuals confirmed that both HexaPro and 2P protein are equally able to detect IgG, IgA, and IgM responses in all sera analysed. Moreover, we extend this observation to include an analysis of glycan engineered S protein, whereby all N-linked glycans were converted to oligomannose-type and conclude that serological activity is not impacted by large scale changes in glycosylation. These observations suggest that variations in glycan processing will not impact the serological assessments currently being performed across the globe.

Mutation-induced changes in the receptor-binding interface of the SARS-CoV-2 Delta variant B.1.617.2 and implications for immune evasion.

Following the initial surges of the Alpha (B.1.1.7) and the Beta (B.1.351) variants, a more infectious Delta variant (B.1.617.2) is now surging, further deepening the health crises caused by the pandemic. The sharp rise in cases attributed to the Delta variant has made it especially disturbing and is a variant of concern. Fortunately, current vaccines offer protection against known variants of concern, including the Delta variant. However, the Delta variant has exhibited some ability to dodge the immune system as it is found that neutralizing antibodies from prior infections or vaccines are less receptive to binding with the Delta spike protein. Here, we investigated the structural changes caused by the mutations in the Delta variant's receptor-binding interface and explored the effects on binding with the ACE2 receptor as well as with neutralizing antibodies. We find that the receptor-binding ß-loop-ß motif adopts an altered but stable conformation causing separation in some of the antibody binding epitopes. Our study shows reduced binding of neutralizing antibodies and provides a possible mechanism for the immune evasion exhibited by the Delta variant.